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We study the plasmon excitations and the electromagnetic response of the magnetically doped 
single layer and multilayer of "helical metals", which emerge at the surfaces of topological insulators. 
For the single layer case, we find a "spin-plasmon" mode with the rotating spin texture due to 
the combination of the spin-momentum locking of "helical metals" and the Hall response from 
magnetization. For the multilayers case, we investigate the electromagnetic response due to the 
plasmon excitations, including the Faraday rotation for the light propagating normal to the helical 
metal layers and an additional optical mode with the frequency within the conventional plasmon 
gap for the light propagating along the helical metal layers. 
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I. INTRODUCTION 

Topological insulators (TIs) [11-11] , as a class of mate- 
rials with insulating bulk states, but conducting surface 
states, have been theoretically proposed [5l-[l3l| and ex- 
perimentally realized both in two dimensions, such as 
HgTe quantum wells 1J|, and three dimensions [la], such 
as Bi2Se3 family of materials [l6l ll7| . For three dimen- 
sional (3D) TIs, for example, Bi2Se3, the surface state 
consisting of a single two dimensional (2D) Dirac cone 
with the gapless Dirac point protected by time reversal 
symmetry since the two components of the Dirac cone are 
Kramers' partners. For simplicity, we can regard these 
two components as "spin" . The 2D surface state exhibits 
a novel "spin-momentum locking" [^ [l^ 11^ 11^ , in the 
sense that the "spin" is locked and points perpendicularly 
to the momentum, forming a left-hand helical texture in 
the momentum space. Therefore, the surface state of TIs 
is also dubbed "helical metals" (HMs)[ll,l23]. For a finite 
Fermi energy Ef, the HM possesses a plasmon excitation 
similar to graphene with the dispersion ^ fc^/^ in the long 
wavelength limit (the momentum fc — > 0). But different 
from graphene, this plasmon excitation is coupled to the 
spin wave due to the helical nature, leading to a novel 
"spin-plasmon" mode(20j. 

Another intriguing feature of a HM is the Hall re- 
sponse due to the lift of the degeneracy at the Dirac 
point by the breaking of time reversal symmetry, when 
the surface of TIs is coated by a ferromagnetic layer, 
or doped with magnetic atoms, which provide magne- 
tization normal to the surface. If the Fermi surface 
lies inside the magneto-bandgap, the Hall conductance 
is half quantized as e^/(2/t), known as "half quantum 
Hall effect" [1, i, i, d iHl • In case of a finite Fermi sur- 
face, the Hall conductance is no longer quantized, but 



still non-zero. For a non-zero Hall response, the longitu- 
dinal and transverse electromagnetic self-excitations are 
coupled, forming a magnetoplasmon, similar to that in 
the conventional metals in a strong magnetic field [2^. It 
is interesting to ask what is the behavior of the "spin- 
plasmon" mode in the magnetically doped HMs. 

In this paper, we investigate the interaction between 
electromagnetic fields and HMs on the surface of a TI 
with magnetic doping, focusing on the plasmon excita- 
tion. In Sec. II, we derive the correlation function for 
a 2D single layer of HMs and discuss how the "spin- 
plasmon" mode is affected by the magnetic doping. The 
electromagnetic waves exist in three dimensions and can 
not interact strongly with the 2D HMs, so we consider 
a model with multilayers of HMs, with an alternating 
stacking of TIs and normal insulators. Based on the di- 
electric function of this model, we discuss the propagat- 
ing modes of the electromagnetic waves in HMs, in Sec. 
III. The conclusion is drawn in Sec. IV. 



II. A SINGLE-LAYER OF HELICAL METAL 

A. Correlation functions 

We start from the calculation of the correlation func- 
tion of a single layer of HMs on the surface(xy plane) 
of TIs with magnetic doping. The single particle 
Hamiltonian [lol l23| for this system is written by 



= hvf{k X (7)2 
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where {cr^;, Uj,, cr^} are the Pauli matrices, ctq is the 2x2 
identity matrix and Vf and Ef are the Fermi velocity 
and Fermi energy respectively. The first term on the 
right hand side is the kinetic term, taking the form of 
Rashba spin-orbit coupling, while the second term gives 
the Zeeman type of coupling due to magnetization. Here 
we only consider the magnetization along z-axis for sim- 
plicity and the Zeeman coupling strength m gives the gap 
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of the HMs. The inplane current operator reads 
7 = evf{z X (j). 



(2) 



According to the Unear response theory, the current- 
current correlation functions read 
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where A is the area of the surface, a, b are x or y, and /? is 
the inverse temperature. Go(fc, is the single particle 
Matsubara's function 
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with the Matsubara frequency a;„. In the long wave- 
length limit fc — >■ 0, the rotation symmetry is restored 
and the correlation function is restricted to the following 
form 



(5) 



where we only keep the momentum independent term 
for both the diagonal part Xo(w) and the Hall response 
Xni'^)-! which are calculated straightforwardly as 
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with X = huj/{2Ef). The step function 0{x) in the imagi- 
nary part of f{x) indicates the appearance of the particle- 
hole excitation for the frequency hu > 2Ef in the fc = 
limit. 

Eqs. ([S]), ^ and ([7]) are the main results of this sec- 
tion. These expressions can recover early results in some 
limits. For example, without magnetization (m = 0), the 
system has a longitudinal plasmon excitation with dis- 
persion proportional to fc^^^ determined by Eq. ([6]), the 
same to the results given in [23|(see next subsection). 
The Hall conductance in the low frequency limit (lu — > 0) 
can be determined by Eq. ([7]) 
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As Ef — > m, the system becomes insulating, and the 
Hall conductance tends to be e^/ (2h), yielding the "half 
quantum Hall effect" jl, i, S [H [2l| . 



B. "Spin-plasmon" modes 

Now let's consider the plasmon excitation in the pres- 
ence of Hall response from Eqs. © and ([7]). Due to 



the existence of the Hall response, the plasmon excita- 
tion is not purely "longitudinal" anymore, but turns out 
to be a mixture of both the longitudinal and transverse 
modes. To calculate these hybrid excitations, we first 
note that for charges and currents distributed in a 2D 
sheet at z = 0, the electric and magnetic fields in the 
plane can be calculated through the Coulomb's law and 
Biot-Savart law 
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where f, f*, J and E are all inplane vectors, and and 
ed are the permeability and permittivity of the surround- 
ings respectively. Here, the magnetic field in the plane 
z = has only z component generated by the inplane 
current. In Eqs. ((TO)) and (fTTj) we neglect the effect of 
retardation 2^ which does not affect the key results dis- 
cussed below and will be considered in the appendix. 

Without loss of generality, we assume the wave is prop- 
agating along y axis in the plane, so that the subscripts x 
and y correspond to the transverse and longitudinal com- 
ponent, respectively. In fact Eq. (fTU]) gives the longitu- 
dinal electric field, while Eq. (fTTj) leads to the transverse 
one according to the Faraday's law (Vxi?)^ = —c~^dtBz- 
Then, after the Fourier transformation of Eq. (fTO]) and 
Eq. (jTI]) . one obtains the inplane electric fields in the 
plane. 



Ex{k,uj) = — — J^(fc,cj) 

C rv 

Ey{k,uj) = Jy[k,uj), 



(12) 
(13) 



where we have applied the continuity condition 
ujp{k,uj) = k ■ J{k,uj), to replace the charge density in 
Eq. (fTU)) by the longitudinal current. Combined with 
Ohm's law J = a{uj)E = i{x/uj)E, Eqs. and ^ 
give rise to the self-consistent equation for electric fields 
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= 0. (14) 



Obviously, without magnetization, there is no transverse 
solution of the above equation, i.e., E^ = 0- Only the 
longitudinal excitation exists with dispersion determined 
by 



27rxo(w) 



(15) 



In the limit lo <^ Ef, one immediately gets the plas- 
mon frequency ui ~ ak^^"^ with the coefficient a = 
\evfkf\/ {^2edEf), recovering the results in [2^. Here 
kf IS the Fermi momentum, and other quantities are de- 
fined as before. It is emphasized that the plasmon disper- 
sion thus obtained has the divergence of group velocity 
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as A; — >■ 0, which is unphysical due to the neglecting of 
the effect of retardation which is rather strong for smaU 
k (see the appendix) . 

In the presence of magnetization, the dispersion of the 
cohective mode in HMs is determined by vanishing the 
determinant of the matrix in Eq. (|14[) , which in the smah 
k limit has the following form 
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In this case both and Ey are nonzero, and the ratio 
between them is determined by 
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It is a pure imaginary number, therefore, the electric field 
is not simply longitudinal but elliptically polarized in the 
plane of HMs for a given k. 

Due to the spin-orbit coupling, the collective mode in 
HMs can also be viewed as the _spin plasmon[20j. Unlike 
the nonmagnetizcd case in Refl20l where a simple inplanc 
linear polarized spin density fluctuation is found, the spin 
density as a vector field in the magnetized case is rotating 
elliptically in the xy plane, and the ratio between the spin 



densities Sx and Sy satisfies 
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The spin texture S{f) in the xy plane is shown in Fig.[T] 
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FIG. 1. An illustration of spin orientation indicated by the ar- 
rows in the xy plane with Sy/ (iSx) positive, and the wavevec- 
tor k is along y axis. 



III. MULTILAYERS OF HELICAL METALS 



A. Dielectric functions and plasmon modes 



FIG. 2. Illustration of layered HMs by stacking topological 
insulator and normal insulator alternatively. In the text, we 
study two basic situations of incidence, either along y or 2: 
axis as indicated by the blue lines. 



Fig. [21 with the magnetization perpendicular to HM lay- 
ers. 

The planes of HMs are equal spacing with the distance 
Id between adjacent layers. The local current density 
can be written as Jioc{k, z,uj) = J„(A;, cj)(5(2; — n/^), 
which is related to the local electric field by the Ohm's 
law Jn{k,uj) = (ix/uj)Eioc{k, z = nld^uj) where the mo- 
mentum k lies in the xy plane. The macroscopic cur- 
rent density and electric field in the x, y direction can be 
obtained by averaging over the local current and local 
electric field, respectively. 
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(19) 



with a test function ^(z — z') following the standard av- 
erage method for macroscopic electrodynamics [l^. We 
assume Id to be much smaller than the wavelength of 
the incident wave, but Id is still large enough so that no 
tunnelings occur between two adjacent layers. 

Under these assumptions and using the definition of 
macroscopic displacement current dtD — dtE + 47r3 we 
obtain the dielectric function 
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where xo and xh are defined in Eq. ([6]) and ([7]) . Eq. ((20|) 
immediately leads to two branches of inplane collective 
modes determined by det[e(a;)] = 0. In the nonmagnetic 
case without Hall response, one obtains the familiar plas- 
mon excitations in the normal metal in the limit lu <ti Ef 
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We now consider the response of the macroscopic elec- 
tromagnetic fields for layered HMs by stacking alternat- 
ing layers of the TIs and normal insulators as shown in 



where n = k'j/{Anld) is the 3D electron density. Since 
electrons can not move in the z direction, the E^ com- 
ponent of electromagnetic waves docs not couple to the 



4 



HMs. Eqs. ((20|) and ([2T|) arc obtained under the assump- 
tion that the electrons in different layers within the range 
of wavelength are moving in phase, since the macroscopic 
fields are the simple average of local ones, which is con- 
sistent with our purpose of discussing optical properties 
for long wavelength of electromagnetic wave. There is 
also another type of excitation in which the electrons in 
the adjacent layers moves out-of-phase 2^ and we will 
not consider it here. In case m 7^ 0, the frequencies of 
the two branches of plasma split and are determined by 



47r[xo(a;p±) ±xg(^P±)] 



(22) 



The electric fields for these two plasmon modes are not 
oscillating linearly, but rotating clockwise and anticlock- 
wise in the plane, respectively. These inplane plasmon 
excitations affect the electromagnetic wave in the bulk 
with inplane polarized electric fields as shown in the next 
section. 



B. Electromagnetic response 



and right-hand modes as seen in Fig. [3l which are ujp^± 
given by Eq. ([22]). In our parameter setting, we have 
Wp-i- > LUp^. Below Wp-, no electromagnetic wave can 
transmit. If ujp- < w < only the right-hand polar- 

ized one is allowed to transmit. If w > the Faraday 
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FIG. 3. Left panel: the dispersion for circularly polarized 
light propagating along z axis, where we take Ef = 0.2eV , 
Id = lOOOA, m = O.SEf and ed = fid = 1. The dispersions for 
left-hand(blue line) and right-hand(green line) polarization 
split. Right panel: the angle per unit length k (rad/mm) is 
plotted as a function of uj. 



Given the wavevector k in normal metal, one can dis- 
tinguish the longitudinal plasmon mode {E \\ k) from 
the transverse optical mode(i? _L k), both of which are 
determined by solving the Maxwell equation 

,2 



{fie)i-j - k 5ij + kikj 



Ej{k,u)=0. (23) 



However in the layered HMs with permittivity given in 
Eq. (PUI) . the Hall response might mix the longitudinal 
and transverse modes, giving rise to the unusual hybrid 
optical modes. Since we are only interested in the influ- 
ence of the collective excitations on the propagating elec- 
tromagnetic wave, we consider two different situations 
with the inplane polarized electric field, one of which is 
propagating along z axis, and the other is propagating in 
the xy plane. For the electromagnetic wave with the out- 
of-planc polarized electric field, it is not affected by the 
inplane collective mode, thus will propagate in a similar 
way as that in a normal insulator. 



1. Propagating along z axis 

In this case we can take k^ = ky = 0, k^ ~ fc, and 
= 0, and the eigcnmodcs of an electromagnetic wave 
arc simply the circularly polarized wave, Ey = ^iE^, 
with the minus sign for the left-hand and the plus sign 
for the right-hand polarization, respectively. Both modes 
arc transverse, i.e., i? _L fc, as in the normal case, but the 
corresponding dispersions are split as 

2 47r[xo(a;±)±Xg(^±)] , k^c^ 
^± — i 1 • 

due to the Hall response, as plotted in Fig. [H There 
are threshhold frequencies at = for the left-hand 



rotation takes place. The Faraday rotation angle can be 
calculated straightforwardly as following [H] 
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(26) 



We plot K as a function of uj in the right panel of Fig. [3l 
where one find n reaches the maximum as w — >■ Wp+ or 
Lj — > 2Ef /h. The giant Faraday rotation at is due to 
the small velocity of the right-hand polarized light at its 
plasma frequency, while the rapid increase of the Fara- 
day rotation angle around hjj/{2Ef) ~ 1 is because the 
real excitation from the lower Dirac cone to the unoc- 
cupied upper Dirac cone above Fermi energy starts to 
happen. The Faraday rotation effect of the magnetized 
surface states of TIs has been carefully investigated for 
the insulating phase [27l[28|. and here we extend the dis- 
cussion to the metallic phase, which is more relevant to 
the present experimental situation [29|. 



2. Propagating parallel to xy plane 

In this case, without loss of generality, we can take 
the incidence along y-axis, i.e., k^ = kz = 0, ky = fc, 
and Ez = 0- The frequencies of the bulk excitations are 
determined by solving the following equation 
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Obviously, if there is no magnetization, i.e., xh = 0, 
the two eigenmodes are decoupled. One of them is the 
longitudinal plasmon excitation, and the other is the 
transverse mode which couples to the electromagnetic 
wave and determines the electromagnetic response of the 
medium. However in the presence of magnetization, the 
Hall response couples these two excitations, leading to 
two hybridized modes in the bulk. The dispersions of 
the two modes are plotted in Fig. IH where we denote the 
upper branch with and the lower one with w_(fc). 

The lower branch, which originates from the longitudi- 
nal plasmon in the nonmagnetic case, has a dispersion 
and affects electromagnetic waves as well due to the mix- 
ing with the transverse mode by the Hall response. In 
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FIG. 4. The dispersion of the two modes of non-transverse 
wave with Ef = 0.2eV, m = 0.5^;/, U = lOOOA, and = 
Md = 1- 

Fig.m one finds two forbidden bands for the transmission 
of the electromagnetic wave, which fall into the regimes 
[0, a;_(0)] and [a;_ (cxj), a;+(0)], where w_((X)) can be iden- 
tified with the plasmon frequency given by Eq. (|2ip 
in the nonmagnetic case, while w±(0) can be identified 
with ijjp± as given in Eq. (|22|) . In the small lo limit, these 
threshold frequencies have the following form 
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FIG. 5. The ratio Ey/{iE^) (left panel) and S'y/(jS'^) (right 
panel) of the inplane elliptic polarization of electric fields of 
the propagating wave with Ef = 0.2eV, Id = lOOOA, m = 
O.SEf and ed = i-id ~ 1- The wavevector k is along the y axis. 
The inset of left panel gives the two inplane elliptical rotation 
modes of the electric fields. 



current accompanied with the electric field E, which can 
be obtained from Ohm's law and Eq. (P7)l 




-E,, 
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Due to the spin-orbit coupling, the ratio between Sy and 
Sx can then be written as 

— = -— = ( (31) 
Sx Jy \edfidi^^ J Ey 



We plot Sy/{iSx) as a function of k in the right panel 
of Fig. [5l which is a real number. This indicates the 
spin orientation in space is rotating circularly at A: = 
0, and cUiptically when fc > 0. This is similar to the 
two dimensional case as shown in section II, the only 
difference is that we have two elliptically polarized wave 
here instead of one. Our calculation shows that one may 
create the spin wave plasmon in stacked HMs by light 
with appropriate frequency and incidence along the xy 
plane. 
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Similar to the single layer case, the electric field is nei- 
ther perpendicular nor parallel to the wave vectors, in- 
stead, it is elliptically polarized in the xy plane, as shown 
in the inset of the left panel of Fig. \5\ where we also 
plot Ey/iiEx) as a function of k which is determined by 
Eq. (HID 



iEx 



(29) 



One can find the wave is circularly polarized at fc = 
where \Ex\ = \Ey\, and as k increases, the wave is el- 
liptically polarized. The upper and lower branches have 
opposite hclicity, and \Ex\ > \Ey\ for upper branch and 



\Ex 



< \Ey \ for lower branch. One can now calculate the 



IV. CONCLUSION 

In conclusion, the plasmon excitation and the electro- 
magnetic response of a single layer and multilayers of 
HMs are carefully investigated for the case with magne- 
tizations. We find that the "spin-plasmon" modes dis- 
cussed in Ref.H^I are modified, with the corresponding 
electric fields, as well as the spin orientation, becom- 
ing elliptical due to the Hall response. Since a single 
layer of HMs can not strongly couple to a 3D electromag- 
netic wave, we consider the electromagnetic response of 
the multilayers of HMs stacked by TIs and normal in- 
sulators. For electromagnetic waves incident normal to 
the HM plane, we find different plasmon frequencies for 
the left-hand and right-hand circularly polarized waves. 
Therefore, for the light with the frequency between these 
two plasmon frequencies, only one type of circularly po- 
larized wave can propagate along the sample. For the 
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frequency above both the plasmon frequencies, a giant 
Faraday rotation is expected. For the hght propagating 
in the HM plane, a new branch of mode appears below 
the conventional plasmon frequency due to the mixing be- 
tween longitudinal and transverse modes. Interestingly, 
the polarization of these two modes are both elliptical 
in the helical metal plane, no longer perpendicular to 
the wavevector k. The new optical modes are expected 
to be observed in an optical transmission and reflection 
experiment [29l - l3l| . 
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Appendix A: The retardation effect in the two 
dimensional plasmon excitation of helical metals 

In this appendix, we consider the retardation effect, 
which has been neglected for simplicity in Eqs. (|12p and 
([T3)) in the text. For this aim, one needs to solve the full 
Maxwell equations subject to the time-varying sources. 
It is then convenient to use the potentials tp and A, which 
satisfy the inhomogeneous wave equations 
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fix, y, z, t) = — p{x, y, t)S{z) (Al) 
A{x,y,z,t) = ^J{x,y,t)SizXA2) 



where c' = c/y^mid, and the Lorenz gauge in dielec- 
tric material V • ^ -I- fid^dC^^dt^ = is adopted. Since 
the system is translation invariant in the xy plane, both 
Eq. (|A1|) and (|A2|) can be written in the following form 
after Fourier transformation 



X{k,z,u)^ f{k,w)S{z), (A3) 



where k is the wavevector in the xy plane, and X and / 
are the potential and source respectively. This equation 
has extended solutions propagating in the full space with 
the general form 



X{k,z,uj) = ae'''^' +[3e- 



■ikz z 



lih^sHkMm 



where k^ = y {lo/c')^ — fc^, and a and /3 are arbitrary 
numbers. There is also a localized solution propagating 
only in the x — y plane which decays exponentially along 
the z direction, 



X{k, z, to) 



f{k, 



-A|z 
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where A = y k"^ — (w/c')^. In this paper we are only 
interested in the localized solution Eq. (|A5p . By substi- 
tuting X with ip and A, and / with p and J, we have 



ip{k,z,uj) ^ e I I, 

cA 



(A6) 



Then the inplane electric ficld(£' = ~Vip~c ^dtA) reads 



E{k,z = 0,cj) 



-p{k,cu) + —^Jik,io) (A7) 



iedX 



Now if one assumes k = (0, k, 0) along the y axis, keeping 
the lowest order terms expanded in terms of c~^, and 
using the continuity condition for the currents, then one 
can obtain Eqs. (Htl) and (fTSj) . 

The group velocity in the nonmagnetized case can be 
calculated to be proportional to fc^/^ as fc — >■ without 
considering the effect of retardation. It is then obvious 
that Vg (X k^^/"^ which is divergent as fc — > 0. This un- 
physical feature can be cured by taking the retardation 
effect, i.e., by using Eq. ()A7p . In the nonmagnetized 
case, the longitudinal and transverse fields are decoupled. 
Eq. (|A7p indicates that the longitudinal electric field is 
determined by the longitudinal current 



(A8) 



According to the Ohm's law Jy = ixoEy/uj, then we 
require 



A = 



27rxo(t^) 



(A9) 



For small w, we can approximate xo(w) ~ a^erf/(27r) with 
a defined in Sec. UbI Then 

—— (AlO) 



a 



4fc2 



If /c 3> a/c', we recover the plasmon frequency in 
Ref. [2^ . As fc — )- we can not neglect the effect of retar- 
dation anymore, in fact, in the opposite limit k <^ a/c', 
we obtain ui sa c'fc. 
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